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Controversy exists regarding the possible existence of a transition between the liquid and glassy 
states of water. Here we use experimental measurements of the entropy, specific heat, and enthalpy 
of both liquid and glassy water to construct thermodynamically-plausible forms of the entropy 
in the difficult-to-probe region between 150 K and 236 K. We assume there is no discontinuity 
in the entropy of the liquid in this temperature range, and use the Adam-Gibbs theory - which 
relates configurational entropy to dynamic behavior - to predict that dynamic quantities such as 
the diffusion constant and the viscosity pass through an inflection where the liquid behavior changes 
from that of an extremely "fragile" liquid to that of a "strong" liquid. 
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It has been evident for at least three decades that 
there is a problem connecting the thermodynamic be- 
havior of water at normal temperatures to that of glassy 
water (which is found below a glass transition tempera- 
ture T g w 136 K) Some contributions assert that 
the liquid and glassy states should be thermodynamically 
continuous, based on measurements of, e.g., specific heat, 
entropy, and relaxation times l^-Q ■ Other contributions 
argue that the liquid and glass should be thermodynam- 
ically distinct j^S-H. In particular, refs. [|||l0| focus on 
a thermodynamically-plausible form for the entropy, and 
determine the limits on the entropy of the glass that are 
consistent with the possibility of continuity. The entropy 
of the glass was subsequently measured |?J , and found to 
be consistent with (but does not require) thermodynamic 
continuity between the liquid and glassy states. 

As a simple illustration of the utility of thermody- 
namics to identify the existence of a transition or other 
anomalous behavior of thermodynamic properties, sup- 
pose we examine the experimental data for the enthalpy 
H, entropy S, and specific heat Cp of liquid water at, 
e.g., 10°C and of ice Ih at, e.g., -10°C. The only way 
to reconcile the large differences in H and S is to hy- 
pothesize a discontinuity in 5 or a large "spike" in Cp 
in this temperature range, even without knowledge that 
a first-order melting transition occurs. In other words 
using only thermodynamic data, one can place relatively 
stringent limits on the thermodynamic behavior near the 
melting transition. Inspired by this fact, we perform test 
for the presence of a "dramatic change" of the liquid ther- 
modynamics below the homogeneous nucleation temper- 
ature Tjj of supercooled water and above the crystalliza- 
tion temperature Tx of glassy water. Specifically, we use 
experimental data on the specific heat, entropy, and en- 



thalpy in both the liquid and glassy states to construct 
two thermodynamically-plausible forms of the entropy in 
the difficult-to-probe region between Tx = 150 K and 
Th = 236 K at 1 atm |llj. We then use these forms for 
the entropy, in conjunction with the theory of Adam and 
Gibbs ]l2fl , to predict behavior of the diffusion constant 
and the viscosity. 

To determine a reasonable form for the entropy 
S = S(T,P), we first focus several of the thermody- 
namic properties that facilitate calculation of S in the 
experimentally-accessible region and also place strict lim- 
its on the possible behavior of S in the difficult-to- 
probe region Tx < T < Th- We define the excess 
enthalpy H cx = H liquid - H clysta i, the excess entropy 
S cx = Siiq U id — S'crystai, the difference of the liquid and 
crystal entropies, and the excess specific heat 
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Each of these quantities is known outside the difficult-to- 
probe region, and in particular at the bounds Tx and Th 
of the experimentally difficult-to-probe region (Table Q). 

• T > Th- H cx (Th) = ffliquid(?ff) — i?crystal(^H) has 

been measured from the heat of crystallization of super- 
cooled water [jl6| 
to S ex by integrating Eq 



We can relate measured values of Cp x 



(0), 



S* CX (T) = 5 CX (T M ) - / -^dT 
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where S cx (Tm) = ASp, the entropy of fusion. We numer- 
ically evaluate the integral in Eq. (||) for T > Th, since 
we know Cp quld from recent bulk sample studies at tem- 
peratures from Tm down to — 29°C [p5|, (and by emulsion 
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studies down to — 37° C and we know Cp ryst for all 

T < T M @. 

• T < Tx- H cx (Tx) = #ii qu id(r\r) — H cl - ysta i(Tx) has 
been measured from the heat of crystallization of glassy 
water |Q . Cp x below Tx is known to be very small, and 
may be taken to be nearly T-independent for T w Tx- 
S e x(Tx) is known from the vapor pressure experiments 
on the glass and and the crystal states Q| . 

• Tx < T < Th- We construct two possible forms 
for for Tx < T < Th similar to the methods 
of refs. §0, but we now include the known value of 
Sex(Tx)- S GX and Cp x fix the endpoints and the slopes 
of S cx at Tx and Th, while the identity 



H ex (T H ) - H ex (T x ) = / C]?dT = 2910 ± 30 J/mol 
Jt x 

(3) 

constrains the area bounded by Cp x - and so also con- 
strains the area bounded by S'cx, due to Eq. ([!]). 

Using these five thermodynamic constraints |l9]| , we 
construct two possible forms of S ox [Fig. |l|(a)] and the 
corresponding Cp x [Fig. [l](b)] . Curve 1 shows the case of 
continuity with no transition, while curve 2 shows conti- 
nuity with a previously-discussed A-transition E 13|. We 
obtain curve 1 in an ad-hoc fashion that satisfies the ther- 
modynamic constraints. We use a closed form for Cp x 
and S'cx in the A-transition case, assuming mean-field be- 
havior 



A more gradual change in S'cx below Th than shown in 
Fig. ^(a) would bound an area J T H S cx dT larger than 
422 J/mol. Furthermore, the inflection in S cx [Fig. [|] 
must occur at T > 215 K, as moving the inflection to a 
significantly lower temperature would also yield an area 
too large. 

We do not hypothesize a discontinuous form for the 
entropy because the magnitude of the possible disconti- 
nuity is unknown. While the data in Table || and the 
thermodynamic constraints do not require a discontinu- 
ity in S con { below Th, the data can not rule out the 
possibility. Furthermore, behavior of Cp x resembling a 
step-function is also possible. However, the accelerating 
increase of Cp x approaching the inaccessible region from 
above makes the A-transition a more natural choice. On 
the other hand, we emphasize that using the available 
data, it is impossible to distinguish the correct from of 
Cf> x , or the exact location of any possible transition (only 
that an anomaly occurs within the approximate range of 
200 K < T < 230 K). 

We now consider the possible implications of the ap- 
proximate forms for S'cx on the dynamic behavior below 
Th- The entropy-based Adam-Gibbs theory |l2| has been 
used to describe the relaxation of liquids approaching 
their glass transitions pl[ |, and provides an explanation 
for the variation of diffusion constant D and viscosity rj 
(or other characteristic dynamic quantity) in the anoma- 
lous range below — 20°C [E2|. We use the prediction 
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from which we obtain 
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In principle, the five free parameters (So, a, b±, and 
T\) may be determined by the five thermodynamic con- 
straints. However, such a procedure yields an experimen- 
tally unreasonable value of T\ = 269 K. To obtain reason- 
able values for the parameters, we choose S ex (200 K) = 
3.4 J/(K-mol), C*p x (200 K) = 15 J/(K-mol), and T\ = 
225 K. This fixes the remaining free parameters in 
Eq. (4) f§. 

Curves 1 and 2 for S cx and Cp x both show sharp 
changes in their behavior just below 230 K. Note that a 
significantly less sharp change in S ex than shown would 
not satisfy the constraint of Eq. (^) , as can also be seen 
by rewriting constraint (||) in terms of the area bounded 
by Sex- From integration by-parts, we find 



Here A is a constant [^3[ . The configurational entropy of 
> t^e liquid, S CO nf = Sii qu id — S v ib , can be understood as the 
^ entropy attributable to the various basins the liquid can 
sample in the energy landscape picture p],p4|. The vi- 
brational component S v ib of the entropy is attributable 
to the thermal excitation the liquid experiences in the 

i i j . ^liquid ^crystal 

basm sampled, so we may approximate S vi ^ « S vi £ 
For typical crystals, S con f ~ 0, since the crystal sam- 
ples a negligible number of basins. Hence S cr y S tai ~ S v ib? 
from which it follows that for the liquid S con f ~ S ex . 



ScxdT = 



TS C 



= 422 ±30 J/mol. (5) 



S^b StaJ is quite good for 



The approximation S^ uld 

liquids near T g |2j|, as the liquid typically samples only 
the deepest basins in the energy landscape, with small 
harmonic excitations about the minimum, similar to the 
crystal behavior. At higher T, the liquid explores a much 
greater region of the landscape, and is no longer localized 
in a single basin for an extended time. Thus at higher T, 
the approximation that S|,'^ ) uld « S^ 8 * 3,1 may not work 
as well @. 

In the case of ice, there exists a residual entropy S rcs at 
zero temperature due to proton disorder p9[ . When we 
subtract S cr ystai from Sii qu id, we also implicitly subtract 
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the contribution arising from S ies , which should also con- 
tribute to the available configurations of the liquid, so we 
include S res explicitly 
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We use this form of 5 C onf to predict the behavior of 
j) and D §(} for T < T H . We select parameters @ in 
Eq. (H) to fit the experimental values of D |32| and 77 |3^] 
[Fig. g). The non-Arrhcnius behavior for T > 230 K 
is typical for a fragile liquid pij. As T — > T g , we find 
ry(r s ) « 10 16 Poise. This is roughly 3 orders of magnitude 
larger than f](T g ) expected from experiments 35 1, but is 
not unreasonable considering that a number of approx- 
imations that were necessary, and further that we have 
extrapolated over 14 orders of magnitude from experi- 
mental data covering only one order of magnitude. Fur- 
thermore, excluding S rcs in Eq. M), results in an absurd 
value r){T g ) w 10 43 Poise, emphasizing the importance of 
including S les as part of S con f f||,|]j ■ 

The dramatic change in entropy around 225 K required 
by the physical constraints is reflected by the "kinks" 
that appear in D and ^atTss225K. In contrast to 
the fragile behavior for T > 220 K, the behavior for 
T < 220 is characteristic of a strong liquid |24|] - i.e. Ar- 
rhenius behavior with an appropriate activation energy. 
We find E D /RT g « 28, E v /RT g « 29, larger than the 
expected activation energy E/RT g w 14 for an "ideal" 
strong liquid |24j], but much less than that of a fragile 
liquid (E/RTg « 80 — 100). Furthermore, adjustment of 
the constants in Eq. (|J) such that r](T g ) — 10 13 Poise, as 
expected from experiments j3f|, yields a value of E v /RT g 
much closer to the "ideal" value for a strong liquid. 
These results support the hypothesis that the fragile be- 
havior of water shows a change to strong behavior for 
T < 220 K. Water is also expected to be a strong liq- 
uid near T g ss 136 K, based on measurements of the 
change in Cp near T g , and the width of the glass transi- 
tion (38) . Such a crossover from fragile to strong behavior 
is not typical of liquids, and merits further experimental 
scrutiny j3{|. 
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TABLE I. Thermodynamic properties of water at 1 atm at 
150 K and 236 K. Here, X cx = Xi iqui( j — A crysta i, the excess 
quantity X of the liquid value relative to the ice Ih value. 
The uncertainties of Sex and JT C x are taken from JD^, which 
contains arguments supporting the reliability of the data. 





T x = 150 K T H = 236 K 


Cf- [J/(K-mol)] 
Sex [J/(K-mol)] 
JTex [J/mol] 




69.2 ± 0.5 [L5-17| 


1.7 ± 1.7 (7 

1380 ± 20 [ 


| 15.2 ±0.1 [13 ; 1€,17] 


1§ 4290 ± 20 |(l|^^ 
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FIG. 1. (a) Possible forms for the excess entropy S ex 
in the experimentally inaccessible region. Curve 1 corre- 
sponds to no transition, while curve 2 shows a A-transition at 
225 K. Any thermodynamically plausible form of S C x (with- 
out a discontinuity) can vary only slightly from these form 
(due to the uncertainty in H cx ). The entropy of fusion 
AS F = 21.8 J/(K-mol) for freezing at 273 K is indicated 
by the arrow, (b) Constant pressure excess specific heat 
Cp X = T(dSex/dT)p for the possible forms of Sox shown in 
(a). Curve 2 has a diverging C P X at the A-transition temper- 
ature. 
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FIG. 2. (a) Diffusivity predicted by Eq. (g). The exper- 
imental data (o) for T > 235K are from 133] . The data for 
T < 160 K (□) are from @ (b) Fit of Sex to viscosity us- 
ing the same procedure. Experimental data (o) are from [ fj3| |. 
Both (a) and (b) show behavior expected for a strong liquid 
for T < 220 K - i.e. Arrhenius behavior with an activation 
energy ~ T 9 /3 (in units of kj/mol) [^4|. The insets show the 
quality of the fit in the region where experimental data are 
available. 
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